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Millions of people are travelling from country to country every day. The major threat is disease spreading that 
causes greater health hazards [9]. Dengue fever is one of the emerging threats now a days throughout the world spread by 
mosquitoes. The mosquito, "Aedes Ageypti" performs the work of a carrier (i.e) the medium for transmitting, for the 
spread of Dengue fever (DF). 

In this paper, Stochastic Ordinary differential equation model (SODE) for DF is proposed for the spread of DF. 
One can understand the underlying processes and develop effective prevention strategies 
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Dengue is a reemerging disease mainly in cities throughout the world. This is found in tropical and subtropical 
region both in urban and semi urban areas. Four distinct but closely related viruses cause this disease [1]. These viruses are 
transmitted to humans by infective female Aedes Ageypti mosquito [6]. 

This fever can go up to 39.5 degree to 41.4 degree Celsius, (i.e) 103.1 - 106.52 degree Fahrenheit[3]. It is 
accompanied by headache and muscle or bone pain. There are also other symptoms like, nausea, vomiting, loss of appetite, 
pain in eyes, and rashes on skin. 

They affect the white blood cells first. Dengue illness is mostly confused with the other viral diseases. This fever 
may vanish and then reoccur along with skin rashes [8]. There is no drug available for the specific treatment of dengue. 
This is also called break bone fever. The symptoms of dengue are evident within 2-7 days of infected mosquito bite [1]. 
There is a need for a mathematical model for dengue in order to protect our society against it. 

MODELING THE EPIDEMIOLOGY 

There are many mathematical models for various diseases. Consider the model Susceptible(S), Infection (I), 
Recovery (R) known as SIR by Kermack-McKendrick [10]. SIR is one of the most basic epidemiological models from 
1927[4]. This is the most widely used model for the spread of disease. In 1957 Macdonald [4] improved the model to a two 
dimensional model with one variable representing human as one variable and the other one representing mosquitoes. 

The extension of this model was proposed by Dietz, Molineaux and Thomas [2]. But this SIR model is not an apt 
model for DF, since contact between the infected and susceptible humans is not the only source for new infections. A 
better dengue model should require tracking of both humans and mosquitoes. Hence, both are included in order to increase 
the accuracy of the model. The model proposed by Ngwa Shu [5] was considered to develop the model for DF. 
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Let [SSj represent the number of susceptible humans at any time t. [^f] represents the number of exposed 
humans at any time t. \jF^ represents the number of infectious humans at any time t. [-^C]^ represents the number of 
recovered humans at any time t. \SS\ denotes the number of susceptible mosquitoes at any time t. \EX\ denotes the 
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number of exposed mosquitoes at any time t. |VfJ denotes the number of infectious mosquitoes at any time t. \TP~^ 
denotes the total number of human population at any time t and \TP\ denotes the total number of mosquitoes population 

L Am 

at any time t. 

Hence [TP] h = [SS] h+ [EX] h+ [lF] h+ [R C ] h and [^^M^M^L 



Let the per capita density dependent death and emigration rate for human be Q. h [jTP^j ■ This is equal to the sum 
of density independent part of the death (and emigration) rate for human j^DZDj and the product of density dependent part 
of the death (and emigration) rate for human j^DDD^ and the total human population at any time 't' \TP~^ ■ 

Hence ^TP\) ^[DID\^DDD] h [TP\ 

The dimension for ^DIDj^ and [^®] is Time' 1 . The dimension for ^DDD^ is Humans' 1 x Time~ l and for 
[ddd] is Mosquito' 1 x Time' 1 

Let the immigration rate for human be fp h . The dimension for p h becomes Humans x Time' 1 . 



Let per capita density dependent death rate for mosquitoes be Q m |jT!Pj J . This is equal to the sum of density 
independent part of the death rate for mosquitoes \DID\ and the product of density dependent part of the death rate for 
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mosquitoes \DDD\ and \TP] the total mosquitoes population at any time 't' 
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HenceQ m ([rP]J = [DID] m + [DDD] m [TP] m 



The probability that the mosquito is infectious is the number of infectious mosquitoes divided by the total number 
of mosquito population. 

Hence we have, r J^ . Let r J™-, = j5 m . 



[TP] ' \TP] 
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Let <P m be the number of times one mosquito would want to bite a human per unit of time, provided the human is 
freely available. The dimension for O m is Time~ l . 

Let O^be the maximum number of mosquito bites a human can have per unit of time. The dimension for O^is 

Time' 1 . 

Let and j^AffiJ represent the number of mosquito bites a human can have per unit of time't' and the 

number of times a mosquito bites humans per unit of time 't' respectively. Hence one can easily interpret that 



[MB]\[TP] h ,[TP 
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and 



[MB] ([TP] [TP] ] = '"i M h The dimensions for \MB] and [MB] is Time~ l . 

I ML Jh L JmJ O m [TP] +® h [ rP \ 

The infection rate of human is the product of the number of mosquito bites a human can have per unit of time 't' 
and the probability of transmission of the viral from mosquito to human P mh . P hm is dimensionless. 



Therefore,^ ^MB\[[TP] h {TP] m )P mh Pm . 

In the similar way, one can obtain the infection rate for mosquitoes. Let this be t m . 

Let [-^Cj^ denote at any time t, the number of recovered humans. T m is the product of the number of times a 

mosquito bites humans \MB] and the summation of the force of infection from the infectious human to the recovered 
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where P, be the probability of transmission of infection from an infectious human to 

nm * J 



a susceptible mosquito given that there is a contact between two occurs and P hm gives the complementary probability of 



P hm . P hm dimensionless, Since, P hm is dimensionless.. 



n^r m = [MBUTP] h {TP]\ 
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Let L be the rate at which the exposed humans move 
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to the infectious state. Then ultimately, the average duration of the latent period for humans is given by -J- . The 

dimension of % h is Time~ l . 

Those humans who are infected only can move to the recovery state. Let the rate of those humans who are moving 
from the infectious state to the recovered state be h . Then the average duration of the infectious period for human is given 

by J- . The dimension of 6, is Time -1 . 

Let L h be the rate of the immune period for humans. Then the average duration of the immune period for humans 

is given by -L- . The dimension of L h is Time~ l . 
L h 

In order to investigate about the transmission of the disease in an area, one has choose the values of ^DIDj^ and 
j^DDDj^ appropriately and that should stabilize the human population. The same interpretation can be given 
corresponding to the mosquito parameters \DID\ and \DDD\ also. Let [fi/T?] and [-6/7?] denote per capita birth 
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rate for humans and mosquitoes respectively. 

In this paper, up to 17 parameters and 9 states are defined. All these parameters and the figure, which is acting as a 
flow diagram, satisfy the following system of Ordinary differential equations, which become SODE model for DF. 

^[^]^^^[[n] +L 4[^]-^ (f) [ 55 i- Q 4t rp ]J 55 i 



' Rc h y e n en - l * [ rc \ - q * {[ tp \ l iF \ -< if i 



dt 



d -[SS m ] = [BIR] m [TP] m -T m (t)[SS] m -Q m {[TP])[SS] n 



d_ 
dt 



dti 

Ex m y T m (t)[ss] m -z m [Ex] m -n m {[TP] m )[Ex] m 

l[/Fj^ m [^] m -^([n)[n ^ 

Analysis of SODE Model 

Dividing each state population by the total population one can scale the state population. Let us denote each 
scaled population by small letters, (i.e) ]^'[ rc ] ft >[ e;t: ] m and^lfj denote at time't', the proportion of exposed 

human, infected human, recovered human, exposed mosquitoes and infected mosquitoes respectively. 
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[rp] m =[^ ]m [rp] m -r[D/D ]m+ [DDD] m [7P] J[rp] ffl (B) 



d_ 

dt 

d_ 
dtV 

It can be shown that DF model (B) has exactly one equilibrium point X jm£s = ^0,0,0,A^,0,0, A^J, with no 
disease in the population. 

By setting, the LHS of (B) equal to zero, J^,[rc]^,[ex;] w ,|jf J all equal to zero and solving for 

\TP] and \TPl one can get, [tp]' = ( W -[Z^H M „ jjjgjjg -^ DDD \i\ and r Jp f = [BIR] -[DID] m 



Now, one of the parameters in epidemiological models is the reproductive number R Q . This R Q can be defined as 

the number of infections that would result from one infectious individual (Either human or mosquito) over the infectious 
period, given that all the other individuals are susceptible. 

This can be defined as, R Q = ^ hm ^ mh . where ^H hm and ^ mh are the number of humans that one mosquito 

infects throughout its infectious lifetime if all humans are susceptible and the number of mosquitoes that one human infects 
through the duration of the infectious period, if all mosquitoes are susceptible respectively. This can be written in 
mathematical notation as, 



4, + [DID] m + [DDD] m [TP] m <p m \TP\ + fa [TP] 



, r ji n . ^ \lDID] m + [DDD] m [TP]']' 
+ <Pi. TP , 



= ^ A0„,P,„„[ TP 1, In +[DID\+[DDD\\TPX\ 

^ + [DlD] h+ [DDD] h [TP]J h [TPl + ^ m [TPl\ h * L ^ 

■ P,„„ + PZ>A (% + lDID] m + [DDD] m [TP]] )"' 



It can also be shown that the disease free equilibrium point x nodjs = ^0,0,0, N h ,0,0, N m j, is locally asymptotically 
stable, if R Q < 1 and unstable if R^ > 1 
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CONCLUSIONS 

In this paper, it is proposed to present a mathematical model for DF. The step by step construction of the model is 
illustrated. The model is analyzed by finding equilibrium point and reproductive ratio. The accuracy of this model can be 
improved by introducing more variables. 
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